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Abstract. A Ginzburg-Landau-like functional is proposed reproducing the main 
low-energy features of various possible high-Tc superconducting mechanisms involving 
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specific heat jump at Tc and the superconducting fluctuations above Tc. Comparison 
with existing data suggests, e.g., that the increase of due to the so-called inter- 
layer tunneling (ILT) mechanism of interlayer kinetic-energy savings is negligible in 
optimally-doped Bi2Sr2CaCu208+a;. 
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1 Introduction 



One of the striking systematics of the superconducting critical temperature of the 
cuprate superconductors (HTSC) is its correlation, within each family of chemically 
similar compounds, with the number N of Cu02 planes per unit cell. [1-3] Well-known 
examples of these families are the Ca-spaced series, with formula ACajv-i (CUO2) n, such 
as the so-called La-, Hg-, TI1-, and Tl2-series (where A is, respectively, La2_a;Sra,02, 
HgBa202, TlBa203_5, and Tl2Ba204_5). In all these Ca-spaced series, the following 
rule is observed to be well obeyed, [1-3] at least for low N values, N <3 [4] : 



where T^=^ and are the critical temperatures of, respectively, the single-layered 
and the A'"-layered compounds (at zero magnetic field and for bulk samples with opti- 
mal hole doping concentration). The origin of relation (1) could in principle be sought 
either in some A^-dependence of the parameters involved in in-plane interactions lead- 
ing to superconductivity, or in the existence of c-direction (i.e., inter-layer) interactions 
contributing, at least in part, to the condensation energy. The latter way of thinking 
is consistent with the proposals made by various authors of different forms of inter- 
layer interactions in HTSC saving energy in the superconducting state. [2, 3, 5] For 
instance, eq. (1) was obtained, for low N, by Leggett [3,4] by considering the screening 
of Coulomb inter layer interactions among carriers. Energy savings in the superconduc- 
ting state occur in this approach mainly in the potential energy of the electrons. [3] 
In what concerns eq. (1) and for low N, Leggett's formalism may be expected to ap- 
ply in essence for a broad variety of superconducting mechanisms responsible for the 
c-direction attractive screening. [3] A notable exception is the so-called interlayer tun- 
neling (ILT) mechanism proposed by Anderson, Chakravarty, and coworkers [2]. In 
the ILT model, savings occur in the kinetic energy due to a deconfinement process of 
the Cooper pairs. This is originated by strong electronic correlations that block the 
coherent interlayer tunneling for single particles, but not for pairs. Strikingly, the ILT 
proposal again leads to eq. (1) in spite of the very different origin of the interlayer 
energy savings. [6] 

An experimental test was proposed by Anderson [7] and Leggett [8] to estimate 
in single-layered (i.e., = 1) HTSC the importance of the ILT mechanism. The 
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test involves the condensation energy Eq, obtained from specific heat data, and the 
c-direction magnetic penetration depth at T = K. Anderson [7] argued that data 
existing for the = 1 compounds of the La-series [9] and Hg-series [10] agreed with the 
ILT prediction. In contrast, subsequent measurements of Ac in the = 1 compounds 
of the La-series, [11] Hg-series, [12] and Tl2-series, [11, 13] concluded that the ILT 
mechanism gives a negligible contribution to Eq. Chakravarty et al. [14] then pointed 
that, in addition to these discrepancies between different measurements, the tests could 
be affected by ambiguities in the obtainment of £"0; this aspect has been recently 
answered in part in [15]. Note that, because of the A^ = 1 limitation, the above 
tests have probed only the energy savings due to the interaction between layers in 
different unit cells. However, this extra-cell interaction may be expected to be the 
less significative one for enhancing Tc\ therefore, the failure of a given superconducting 
interlayer mechanism to account for the condensation energy of the A^ = 1 cuprates 
does not rule it out completely as a substantial source for the enhancement of when 
A^ > 1 (see also below). 

In the present work, we introduce a Ginzburg-Landau (GL)-type functional that 
reproduces in essence the main low-energy features of the above proposals of interlayer 
kinetic and/or potential energy savings in HTSC. This functional is based on a simple 
energy-balance argument expected to be a good approximation near the transition and 
for low Cooper-pair densities, and may be useful for finding experimentally testable 
relationships involving the interlayer superconducting energy savings in HTSC with 
A" > 1. Two examples of such tests will be given in this letter: The first is a rela- 
tionship between A" and the mean-field specific heat jump at Tc- If it is fulfilled in 
a family of chemically similar HTSC, it will indicate that the superconductors in the 
series differ only in their interlayer interactions (i.e., that eq. (1) is actually due to a 
superconducting interlayer mechanism). The second test will be provided by the Gaus- 
sian superconducting fiuctuations above Tc- Our present functional leads, for zero or 
weak magnetic fields, to fiuctuations identical to the well-known ones of multilayered 
superconductors with no interlayer energy savings; however, the involved parameters 
acquire now additional meanings. In particular, the interlayer Josephson tunnelings are 
related to the maximum increase of Tc that could be attributed to the ILT mechanism. 
Finally, we end this letter with a brief discussion of available experimental data. 
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2 A preliminary energy-balance argument leading 
to eq. (1) 



In a layered superconductor with N superconducting layers in each c-direction unit cell 
of length s, we may write the free energy A/ saved in the superconducting state, per 
Cooper pair and unit volume, as the following sum: 

A/ = [TVA/ii + {N- 1) A/^int + A/^ext] . (2) 

Here A/y, A/^int and A/^ext ^-re the energy savings, per Cooper pair and unit area, 
due to, respectively, the in-plane interactions in each superconducting layer, the inter- 
actions between each two adjacent layers in the same cell (intra-cell interaction), and 
the interactions between each two adjacent cells (extra-cell interaction). We assumed 
in eq. (2) that all the superconducting layers and intra-cell separations between them 
are equivalent. We also assumed that the energies saved by each inter layer interaction 
may be considered as independent, and that the interactions between non-adjacent lay- 
ers are negligible. These assumptions are expected to be good approximations in both 
Leggett's and ILT proposals. [2, 3] Let us now expand A/y in powers of the reduced 
temperature s\\{T) = la{T/Tc\\) ~ (T — Tc\\)/Tc\\, where Tc\\ is the critical temperature 
of the system if the interlayer interactions were absent; this expansion is obviously 
reminiscent of the type of reasonings used in the GL-like theories. So, for the in-plane 
interactions above S\\{T) = we write A/y = —ae\\{T), where a is constant and posi- 
tive. For the interlayer interactions, as they will be less dependent on £11(7"), we may 
write in first approximation A/^int = OiSi^t and A/^ext = ce^ext, where 5int and Sext are 
dimensionless constants (positive, if the corresponding interactions help superconduc- 
tivity). After these simple power expansions, the actual critical temperature can 
be easily calculated by just writing the condition Af(T^) — 0, i.e., 

£fs(ri^) = o, with £,^ff(r)^£||(r)-^-(i-l)5,.t. (3) 

which directly leads to: 

C C ' ' ^ r- S I ^ ^ \ 



-<N=1 



T 

c c 



In ^ = (5int - <^ext) (l - ■ (4) 



We find, therefore, a result equivalent to eq. (1), independently of whether the mech- 
anism of interlayer energy savings imply the kinetic or potential energies, or both. In 
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the above equation we have considered it useful to exphcitly emphasize the fact that 
the logarithm of the quotient of two temperatures is approximately equal to the rel- 
ative distance between them; we will omit the explicit emphasis of this point in the 
remainder of this paper. 

Note that ior N — 1,2 the above equations also lead to: 

These relationships indicate that if 5cxt "C (^mt is possible for T^='^ to be quite different 
from T^^^ even if T^=^ ~ Tc\\. In this respect, we note that the recent observation 
of superconductivity in a sample of the = 1 La-series compound with thickness 
only one c-direction unit cell [16] is not contradictory with the existence of interlayer 
superconducting energy savings in HTSC. These aspects also enhance the interest of 
testing the existence or not of interlayer superconducting energy savings in the N > 1 
compounds. 



3 A simple Ginzburg-Landau (GL) free-energy func- 
tional for HTSC with interlayer energy savings 

Our next step is to introduce a GL-type functional consistent with the main proposals 
for interlayer kinetic and/or potential energy savings in HTSC. This functional must 
be chosen so that for low energies and Cooper-pair densities it recovers the energy 
balances of the previous section. This is fulfilled by the following expression for the 
difference between the superconducting and normal-state free energies, AFf^^], at zero 
magnetic field: 

r ^ ( h 

AF[*] = /d^rX:Eao|£||(r)|*,g^ + 2^l*,„r + d(o)|v.,*,„r 

-7. {^jn^U,n + c-c.) + (27. - S,) (^l2^l±^l±lfl^ I . (6) 

In the above expression, r = (x, y) and S/xy are the in-plane coordinates and gra- 
dient (we neglect the possible in-plane anisotropy); the indexes {j,n) label each j'th 
superconducting plane of the nth c-direction unit cell (we use also {N + l,n) for the 
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(l,n+ 1) plane); \l/j„ = \l/j>),(a;, ?/) is the superconducting order parameter of the {j,n)- 
plane; ^,abiO) is the GL ampUtude of the in-plane coherence length; and b are the 
GL constants. The interlayer interactions in the above functional are contained in the 
last two summands. In them, 7^ is the usual Josephson coupling constant between the 
(j, n) and (j + l,n) planes, and 5j parametrizes the interlayer energy savings due to 
interactions also between these two planes. These 6j parameters coincide with the ones 
already used in eqs. (3) - (4), and to our knowledge this is the first time that they 
are introduced in a GL model of a layered superconductor. As we did in the previous 
section, we may assume in the HTSC only two different interlayer separations (intra- 
and extra-cell) and so we take jn = 7ext, = ^ext and jj^tN = Jint, ^jj^N = ^mt- 

The parameters 5j deserve further discussion. When they are zero, the interlayer 
summands in eq. (6) reduce to 7j|^jn — ^j+i,nP- As the latter expression is never 
negative, it does not lead to energy savings. In fact, this is precisely the well-known ex- 
pression for the interlayer interaction in the extended-Lawrence-Doniach GL-functional 
for multilayered superconductors without interlayer energy savings. [17-19] In contrast, 
when 5j > the interlayer energy can be negative: For instance, if 5j — 27^ the last 
interlayer summand in eq. (6) is zero and the first one becomes proportional to cos (pj, 
where cpj is the difference of the phases of and ^j+i,n- This is precisely the form 
of the interlayer kinetic-energy savings proposed on microscopical grounds by the ILT 
model. [2,7] Finally, interlayer potential energy savings as those in Leggett's proposals 
for the superconducting mechanism [3] may also be crudely included in our functional 
by considering different Sj > 0. This is coherent with the fact that these savings may 
be expected to arise in the second of the interlayer interaction summands, i.e., the one 
not involving the phase of the superconducting wave function. Let us also note that 
these kinetic and potential interlayer energy savings could coexist in the HTSC; in that 
case, the Sj would result from the sum of the contribution of each energy-saving source, 
leading then to energy savings essentially of the form A + B cos (fj, as in fact proposed 
in [8] for that mixed scenario. 

It is not difficult to calculate the resulting from eq. (6). For that, one just has 
to write the condition AF[\E'j„(0)](Tj^) = 0, where \E'j„(0) are the equilibrium wave 
functions, verifying dAF/d'i/jn — 0. Just as expected, simply follows eqs. (3) - (5). 
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4 Jump at Tc of the mean-field specific heat 



It is also quite direct to calculate from eq. (6) the mean-field specific heat jump at the 
transition, Cj^mp, i.e., the discontinuity at T = in the specific heat when considering 
the uniform wave function minimizing AF[^]. We get: 

- - ^ (7) 



N 

jump 



^jump — ^rpN ^ ' 

where again s is the layered-structure repetition distance. This proportionality of c 
and N/ sT^ has been to the best of our knowledge unnoticed up to now even for layered 
superconductors with no interlayer energy savings. Two important additional remarks 
are: First, cj^mp found to be independent of the interlayer interaction parameters, 
both 5j and 7j (except indirectly trough Tj^, see eq. (4)). Second, according to eq. (7) 
the quantity c^um^/ {N / sTc) will be the same for all the HTSC with equal in-plane 
interactions (and hence equal GL parameters oq and b). In particular, experimental 
verification of this constancy in a series of chemically similar HTSC would provide a 
compelling argument favoring that eq. (1) is in fact due only to the existence of a su- 
perconducting interlayer mechanism, and not to some A^-dependence of the parameters 
involved in in-plane interactions. (We note however that failure of fulfilment of that 
proportionality would not rule out completely interlayer interactions contributing at 
least in part to the Tc enhancement, as they could still coexist with in-plane effects 
varying oq and h in the series.) 



5 Gaussian-Ginzburg-Landau (GGL) fluctuations 

above Tc 

Let us now consider the superconducting fluctuations above Tc that result from func- 
tional (6). For that, we apply to it the same formalism previously used in [18, 19] to 
study the Gaussian regime of the superconducting fluctuations in multilayered super- 
conductors without interlayer energy savings, i.e., with 5j = 0. The calculations go in 
parallel to those in [18, 19], and so we will not make them explicit in this letter. The 
important result is that the fluctuation spectrum of functional (6) with arbitrary 6j 
is just the same as found in [18, 19] for Sj — 0, if we interpret in these equations Tc 
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as . This applies also to the results directly derived from that spectrum, including 
the fluctuation-induced observables calculated in [18, 19] at T > Tc for zero magnetic 
flelds and also for weak magnetic flelds perpendicular to the layers, i.e., the in-plane 
paraconductivity [18], the in-plane magnetoconductivity [18], the fluctuation suscepti- 
bility [18], and the fluctuation speciflc heat [19]. Therefore, for all these observables we 
conclude that the final expressions obtained in refs. [18, 19] (i.e., eqs.(3.4)-(3.8) and 
(4.5)-(4.7) in ref. [18], and (9)-(12) in ref. [19]) are applicable for all values ofSj. 

The above result indicates that any analysis of experimental data based on the 
above equations of refs. [18,19] will be also apphcable in the 5j ^ scenario. This is an 
important information because those analyses are able to determine rather unambigu- 
ously the values of the Josephson couphngs between layers, 7int and 7ext (see [18-22] 
and also below) . Precisely these Josephson couplings give a maximum limit for the in- 
terlayer kinetic energy savings in eq. (6). In particular, the maximum relative increase 
of Tc that may be due to the ILT mechanism is given by eq. (5) with 5int = 27int and 
(5ext = 27ext- Note that in [18-22] instead of the variables 7int and 7ext it is often used 
the equivalent set composed by ^c(O) and 7int/7ext, where ^c(O) is the c-direction GL 
coherence length amplitude; in terms of these alternative variables, eq. (5) leads to a 
maximum increase of due to the ILT mechanism of: 



6 A brief comparison with existing experimental 
information 

Testing eq. (7) in a given family of HTSC requires reliable knowledge of Cjump on it. 
Unfortunately, the latter proves difficult. An illustrative example of these difficulties 
happens in the Tl2-series. The N = 2 and 3 compounds of this series were measured 
by the same group, in similar samples, and using the same criterion for obtaining cjump, 
in [23]. The corresponding c-^ump/ / sTf]) are equal to 2.6 x 10^"^ J/m^ for N = 2, and 
2.8 X 10-3 J/m^ for N = 3. This suggests fulfillment of eq. (7), i.e., that the T^iN) 




(8) 




8 



variation in this series is due only to interlayer interactions. However, measurements by 
other group [24] in the = 1 compound of the same series revealed a quite small and 
symmetric specific heat peak around T^. This may be interpreted as a negligible Cjump 
(and then not fulfillment of eq. (7)), or could have its origin in, e.g., Tc-inhomogeneities 
broadening and symmetrizing the specific heat peak around Tc [25], or in pseudogap- 
induced effects [24]. Therefore, an ultimate testing of eq. (7) would imply a more 
extensive analysis of the experimental specific heat. Naturally, it should also include 
other Ca-spaced series. Such a detailed study is beyond the scope of the present letter. 

Much less ambiguity exists at present in understanding the superconducting fluctu- 
ations above Tc in HTSC in terms of the GGL model of multilayered superconductors. 
As commented above, for our present purposes the main interest lies in the N > 1 
compounds. Measurements in high-quality single-crystal samples are available for 
at least two bilayered and optimally doped HTSC, Bi2Sr2CaCu208+a; (Bi-2212) and 
YBa2Cu307_5 (Y-123). Those data were extensively analyzed in terms of the GGL 
model of multilayered superconductors with null interlayer energy savings. [18-22] As 
mentioned above, these kinds of analyses remain fully applicable in the case of nonzero 
interlayer energy savings, with the bonus that the obtained Josephson coupling param- 
eters give the maximum increase of Tc that could be attributed to the ILT mechanism. 
For instance, in [18,19,22] it was determined for Y-123 that ^c(O) ~ 1.1 ± 0.1 A, and 
1/30 ^ 7int/7ext ~ 30. By using eq. (9), these values suggest an upper limit of around 
70% for the increase of Tc due to kinetic-energy savings. In the case of Bi-2212, the 
results are more conclusive: The analysis of the superconducting fluctuations in this 
compound leads to ^c(O) ~ 0.5 A, and to 7int/7ext values whose boundaries compatible 
with experiments depend on ^c(O) roughly as ^^(0) ~ 7int/7ext ~ Cj^^(O) (where ^c(O) is 
in A). [21,22] When these values are used in eq. (9), they lead to around 1% maximum 
increase of Tc due to the ILT mechanism in Bi-2212. 



7 Conclusions 



We introduced a simple Ginzburg-Landau-type functional that reproduces the main 
low-energy features of the existing proposals of interlayer kinetic and/or potential en- 
ergy savings in HTSC. This functional may be used to flnd relationships between these 
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savings and experimental observables. Two examples of such relationships were pro- 
posed: The first involves the mean-field specific heat jump at Tc, Cjump, and the normal- 
ized mean interlayer distance sTc/N; if both are inversely proportional to each other 
in a series of HTSC, this would indicate that they share identical in-plane interactions 
and their Tc's are different due only to interlayer interactions. The second relationship 
involves the superconducting fluctuations above Tc at zero or weak magnetic fields, 
from which it may be obtained a maximum limit for the relative increase of Tc due to 
interlayer kinetic energy savings (as those in the ILT model). When compared with 
available experiments, this second relationship indicates that the increase of Tc due to 
interlayer kinetic energy savings is negligible in optimally-doped Bi-2212, and between 
zero and ~ 70% in optimally-doped Y-f23. Although based on a simplified model, 
these conclusions can be expected to be, at least, qualitatively correct. They provide, 
to the best of our knowledge, the first test of the significance of the interlayer kinetic 
energy savings in the T^ of any N > 1 HTSC. 



Acknowledgements 

The main part of this research was performed in the Department of Physics of the Uni- 
versity of Illinois at Urbana-Champaign thanks to support from the Fulbright Founda- 
tion and to the hospitality of Prof. A.J. Leggett, who first called my attention to the 
problem of the interlayer condensation energy in cuprates and with whom I enjoyed 
various valuable talks on the subject. I am also grateful to Prof. F. Vidal for many 
enlightening discussions about the superconducting fluctuations and the jump at Tc of 
the specific heat. Financial support is also acknowledged to the CICYT, Spain, under 
grant No. MAT2001-3053. 



10 



References 



[1] Cava R.J., J. Am. Ceram. Soc. 83, 5 (2000). 

[2] Anderson P.W., The Theory of Superconductivity in High-Tc Cuprates (Prince- 
ton University Press, Princeton, NJ) 1997, and references therein. 

[3] Leggett a. J., J. Phys. Chem. Solids 59, 1729 (1998); Phys. Rev. Lett. 83, 
392 (1999). 

[4] The difficulties that appear when trying to apply eq. (1) for > 4 are commented 
in [3] and also in Jansen L. and Block R., Phys. Rev. Lett. 85, 3983 (2000); 
Leggett A.J., Phys. Rev. Lett. 85, 3984 (2000). See also [26]. 

[5] For some examples of other proposals of superconducting mechanisms involving 
interlayer interactions in cuprates see, e.g., Tesanovic Z., Phys. Rev. B 36, 
2364 (1987); Harshman D.R. and Mills A. P., Jr., Phys. Rev. B 45, 10684 

(1992) ; Kettemann S. and Efetov K.B., Phys. Rev. B 46, 8515 (1992); 
LlEGHTENSTEiN A.I., Mazin I.I. and Andersen O.K., Phys. Rev. Lett. 74, 
2303 (1995); SuLE P., cond-mat/0303585. 

[6] Another possible way to obtain eq. (1), at least approximately, is using Kosterlitz- 
Thouless-transition ideas. See, e.g., Matsuda Y. et al, Phys. Rev. B 48, 10498 

(1993) ; Corson J. et al, Nature 398, 221 (1999). 

[7] Anderson P.W., Science 268, 1154 (1995); 279 (1998) 1196. 
[8] Leggett A. J., Science 274, 587 (1996). 

[9] UCHIDA S., Tamasaku K. and Tajima S., Phys. Rev. B 53, 14558 (1996). 
[10] Panagopoulos C. et ai, Phys. Rev. Lett. 79, 2320 (1997). 
[11] Schutzmann J. et al, Phys. Rev. B 55, 11118 (1997). 
[12] KiRTLEY J.R. et al, Phys. Rev. Lett. 81, 2140 (1998). 

[13] MOLER K.A. et al. Science 279, 1193 (1998); Tsvetkov A. A. et al. Nature 
395, 360 (1998). 

[14] Chakravarty S., Kee H.-Y. and Abrahams E., Phys. Rev. Lett. 82, 2366 
(1999); Phys. Rev. B 67, 100504(R) (2003). 

[15] VAN DER Marel D. et al, Phys. Rev. B 66, 140501 (2002). 

11 



[16] Bozovic I. et al, Nature 422, 873 (2003). 
[17] Klemm R.A., Phys. Rev. B 41, 2073 (1990). 

[18] Ramallo M.V., POMAR A. and ViDAL F., Phys. Rev. B 54, 4341 (1996), and 
references therein. 

[19] Ramallo M.V. and Vidal F., Phys. Rev. B 59, 4475 (1999). 

[20] In [18,19] a BCS value was imphcitly assumed for the relaxation time of supercon- 
ducting fluctuations. This was shown to be in fact correct for HTSC in RAMALLO 
M.V. et ai, Europhys. Lett. 48, 79 (1999). 

[21] POMAR A. et al, Phys. Rev. B 54, 7470 (1996), and references therein. 

[22] A review of the analyses of superconducting fluctuations in HTSC in terms of the 
multilayered GGL model without interlayer energy savings was presented in ViDAL 
F. and Ramallo M.V., The gap symmetry and fluctuations in high temperature 
superconductors, edited by J. Bok et al. (Plenum, London) 1998, p. 477. These 
analyses have been recently extended also to the high-£ region {e > 0.1); see, e.g., 
Vidal F. et al, Europhys. Lett. 59, 754 (2002) and references therein. 

[23] JUNOD A. et al, Physica C 159, 215 (1989). 

[24] LORAM J.W. et al, Physica C 235, 134 (1994). 

[25] Ramallo M.V., Carballeira C. and Vidal F., Physica C 341-348, 173 
(2000). 

[26] After the flrst version of the present e-print appeared, it was posted the work 
Chakravarty S., Kee H.-Y. and Voelker K., cond-mat/0309209. There, 
flrst these authors write a functional similar to our eq. (6) with Sj = 2jj (corre- 
sponding therefore to the ILT limit) and then they propose, in order to explain 
the dependence of Tc on hole doping, and also on for > 4 (when eq. (1) may 
fail [4]), to add to that functional the possible effects of a competing (pseudogap) 
order parameter and of charge (doping) imbalances between nonequivalent layers. 
These effects are not expected to signiflcantly influence the tests discussed in the 
present letter, at least for the N <3 and optimally-doped case on which we have 
focused. 



12 



